Absence of supersensitivity to small input signals in generalized on-off systems 
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It has recently been shown that nonlinear skew product dy- 
namical systems with invariant subspaces which are capable 
of displaying on-off intermittency can show supersensitivity 
to small input signals. 

Here we show that this supersensitivity is absent for more 
general dynamical systems with non-skew product structure, 
capable of displaying a generalized form of on-off intermit- 
tency, and is therefore in this sense fragile. This absence of 
supersensitivity is of importance in view of the fact that dy- 
namical systems are generically expected to be of non-skew 
product nature. 



Many dynamical systems of interest possess symme- 
tries or constraints which force the presence of invariant 
subspaces. A great deal of effort has recently gone into 
the study of such systems (sec e.g. A sub-class of 

these models, namely those with skew product structure 
(and normal parameters^]), have been shown to be capa- 
ble of producing a number of novel modes of behavior, 
including on-off intermittency and bubbling |jj . 

Recently, Zhou and Lai Q] have shown that systems of 
this type can display supersensitivity, in the sense that 
small constant or time varying inputs to the system can 
induce extremely large responses. The authors further 
claim that with an additional odd symmetry condition, 
this supersensitivity is robust to addition of noise. Such 
supersensitivity could be of importance in many fields, 
including the study of synchronization of coupled chaotic 
systems || and the design of sensor devices jgj . 

The results on on-off intermittent systems reported 
by these authors can all be described within the frame- 
work of skew product systems. Generically, however, one 
would expect typical dynamical systems to have non- 
skew product structure (with non-normal parameters). 
Systems of this type have recently been studied and have 
been found to be capable of displaying a number of addi- 
tional novel dynamical modes of behavior, absent in skew 
product systems, including a generalization of on-off in- 
termittency, referred to as in-out intermittency M . 
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1 Parameters that leave the dynamics on the invariant mani- 
fold unchanged are called normal, otherwise they are referred 
to as non-normal. 



The easiest way to characterise in-out intermittency 
is by contrasting it with on-off intermittency. Briefly, 
let Mi be the invariant subspace and A the attractor 
which exhibits either on-off or in-out intermittency. If 
the intersection Aq = A n Mj is a minimal attractor, 
then we have on-off intermittency, whereas if Aq is not 
a minimal attractor, then we have in-out intermittency. 
In the latter case there can be different invariant sets in 
Aq associated with attraction and repulsion transverse 
to Aq, hence the name in-out. Another crucial difference 
between the two is that, as opposed to on-off intermit- 
tency, in the case of in-out intermittency the minimal 
attractors in the invariant subspaces do not necessarily 
need to be chaotic and hence the trajectories can (and 
often do) shadow a periodic orbit in the 'out' phases 0. 

Our aim here is to find out whether this type of super- 
sensitivity, observed in on-off intermittent systems, per- 
sists in more general non-skew product systems which 
are capable of displaying in-out intermittent behaviour. 

A simple class of maps that can model both on-off and 
in-out types of intermittency is given by 



_i = F(x n ,y n ,a), y n +i= G(x n ,y n , a), 



(1) 



where G(x n ,0,a) = 0, the variables x n and y n represent 
the dynamics within the invariant submanifold (y = 0) 
and the transverse distance to it respectively and a € TZ m 
are the control parameters of the system. A special sub- 
set of these systems, with skew product form over the 
dynamics in x, can be written as 

x n+1 = F(x n ,a), y n+1 = G(x n ,y n ,a). (2) 

By considering a skew product system of type (^), Zhou 
and Lai [Q] modelled the motion near the invariant sub- 
manifold y — 0, using a Fokker Planck equation. In this 
way they were able to predict that the sensitivity S of 
the map in the neighbourhood of a blowout bifurcation, 
leading to on-off intermittency, is given by 
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(3) 



pln(r/p)' 

where (y) is the average of the transverse variable y, p 
is the input signal and r is the threshold below which y 
goes through a laminar phase. They were able to confirm 
this prediction numerically. 

To study whether non-skew product (in-out intermit- 
tent) systems can also display supersensitivity, we con- 
sidered a particular example of the map (0) in the form 



Xn+1 — ^*^n(l X n ) d" SX n y n: 
bx "-i 

y n+ i = ve n y n + ay n , 



(4) 
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where r £ [0, 4] and (s, v, a,b) € 1Z are the control param- 
eters. Note that for s = 0, the map (|]) has the skew prod- 
uct form ^ and for fixed r, the parameters s, a^boxv are 
normal. Thus depending upon the choice of its param- 
eters, this map is capable of displaying both on-off and 
in-out types of intermittency, Note also that this map 
possesses the odd symmetry condition G(—y) = —G(y) 
that was found in |4[] to be required for the robustness 
of supersensitivity with respect to noise. Also the trans- 
verse Lyapunov exponent At for this map can be readily 
calculated to be 



At = In ; 



(5) 



where (x) r is the average of the variable x„ for an initial 
condition on the invariant submanifold y = 0. 

To study the effects of an input signal on in-out sys- 
tems, we considered a variant of this map given by 



Vn+1 



rx n (l - 
ue hXn y n 



+ ayl + P, 



(6) 



where the real parameter p models the effects of a small 
input signal. 
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FIG. 1. Comparative study of bursting behaviour in 
in-out (left panels) and on-off (right panels) regimes, for an 
input signal p = 10~ 10 . Note how in-out dynamics (lower left 
panel) is insensitive to the input signal, in comparison with 
on-off (lower right panel) . The upper left panel and the inset 
on the lower left panel also demonstrate clearly the presence of 
the period 12 attractor in the invariant submanifold. The pa- 
rameters values are r = 3.8800045, v = 1.82, b = -1.020625, 
a = — 1 and s = —0.3 for the in-out case and r — 3.82786, 
v — 1.82, b = —1.006, a = —1 and s = for the on-off case. 

To begin with, we made a comparative numerical study 
of the sensitivity of in-out and on-off systems to input 
signals p, using @. Fig. |l| shows a comparative study 
of the bursting behaviour in the two cases close to, but 
below, their blowout points. As can be seen from the 
comparison of the lower panels, there is very little burst- 
ing in the in-out case. 



To further demonstrate this relative insensitivity in the 
in-out case, we made a study of the sensitivity S of the 
systems close to their blowout points, as a function of 
the input signal p. This is shown in Fig. 0, which again 
demonstrates a distinct absence of supersensitivity for 
the in-out case, specially for the lower input signal levels. 
Furthermore, it shows a saturation in sensitivity in the 
in-out case for input signals p ^ 10 ~ 7 . 
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FIG. 2. Dependence of the sensitivity S on the input 
signal p for in-out and on-off cases. Note the relative insen- 
sitivity of the in-out to the input signal and the saturation in 
sensitivity in this case. The parameter values are as in Fig. 



These results indicate a clear lack of supersensitivity 
in the in-out case to small constant input signals. 

To understand this qualitative difference between the 
on-off and in-out cases, we briefly recall a number of 
differences between the two cases, relevant to our dis- 
cussion here. In the case of on-off, the attraction and 
the ejection of the orbits near the invariant manifold are 
brought about by a single chaotic attractor in Mj. Thus 
for the values of the control parameter close to but be- 
low the blowout point, the chaotic attractor in Mj be- 
comes transversally weakly attracting, but there can be 
repelling orbits within this attractor that are transver- 
sally unstable, leading to bubbling and allowing the orbit 
to access the lower and upper boundaries frequently. The 
system thus becomes sensitive to small inputs, producing 
large bursts and hence supersensitivity. 

For the in-out case, on the other hand, the 'in' and 
'out' phases are governed by two separate invariant sets 
in Mj: a chaotic saddle and a periodic attractor respec- 
tively. Thus for the values of control parameter (b in 
our case) above the blowout value, the chaotic saddle 
in the invariant submanifold is transversally attracting 
whereas the periodic attractor in Mi is transversally un- 
stable with a positive transverse Lyapunov exponent At- 
As a result, an orbit drawn towards the invariant sub- 
manifold by the chaotic saddle is thus ejected by the 
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transversally unstable periodic attractor, leading to in- 
out intermittency. On the other hand, for the values of 
control parameter b just below the blowout value, the 
unique periodic attractor in the invariant submanifold 
becomes transversally stable (with At < 0), while the 
chaotic saddle still remains transversally attracting. As 
a result orbits drawn towards the invariant submanifold 
by the chaotic saddle get attracted to the periodic orbit 
there (see Fig. || for a schematic depiction) . 



y 2 = f v 2 e Hxo+Xl) + ve bxi + l) p, 



J/12 = I v ™ e Kxo+xi+...+xu) 



. + ue bx ^ + l) p, 

( 
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FIG. 3. Schematic diagram showing the main dynamical 
features of the in-out process near the blowout point with 
At < 0. This picture does not qualitatively change in presence 
of a input signal p, even though the period 12 attractor and 
the chaotic saddle are slightly shifted from the previous Mi 
(represented by the dotted line). 

We shall now show that the presence of a small input 
signal p leaves the above dynamical picture essentially 
unaltered, apart from displacing the location of the pe- 
riodic attractor off the previous invariant submanifold. 
There are two ways to see this. Firstly, for small values of 
the input signal p <C C(l), the periodic orbit is expected 
to persist by continuity. We numerically confirmed that 
the period 12 orbit involved in the in-out intermittency 
studied here (see for details) does indeed survive for 
small values of p, albeit shifted slightly off the invariant 
submanifold Mj (see the left panels of Fig. [[]) . 

Alternatively, we can estimate the transverse location 
of the displaced periodic orbit. To do this we recall 
that we are interested in small displacements from Mi, 
which implies that as a first approximation we may ignore 
higher order dependence on y. We therefore approximate 
the map (ph by 



Xn+l 
Vn+l 



ve bXn y n +p, 



sx n y n , 



(7) 



where the second order term in y has been kept in the 
x map in order to ensure the essential overall non-skew 
product structure of the system. 

The period 12 attractor involved in this case has x n 
values satisfying x n +\2 — x n and y n values given by 

2/o = P, 

2/i = (ve bx °+l)p, 




b\ g [ 1 + 11 - { l™° dl2) \x 



(i + n-l) m odl2 



p, 



where |_ J denotes the integer part. 

The above expressions for y n change periodically (with 
period 12), depending on the initial x. The asymptotic 
average value of y can then be approximated by 



lim (y) w I v j e \-£o 12 j i , K 

ro^oo — ' \ 
j=0 



12 S Xi 

1 — ve i=0 1 p 



1 - ve XT 



At — In 



(8) 



where we have used (||). Interestingly this enables us to 
find the sensitivity S as a function of At 
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(9) 



which is independent of p, thus explaining the saturation 
in sensitivity S observed in the in-out case in Fig. ||, 
in clear contrast to expression for the on-off sensitivity 
given by (|) . 

It now remains to show that apart from the above shift 
off the submanifold, the periodic attractor remains es- 
sentially intact. To see this, recall that the effect of a 
non-zero p on x is given by 



x n+1 ~ F n (xi, a) + sx n {yY 



(10) 



where F n (xi 1 &) represents the x component of the n th 
iterate of the map (0). Using (0) this gives 



\x n {yf ~\s\0{\) 



(11) 



Now for input signals p ^ 0(1) and for the parameter 
regimes chosen here, 1 — e Ar 3> *Jp, which implies 



\s\x n {yf «P, 



(12) 



showing that to this approximation the p-induced vari- 
ations in x are extremely small (relative to p), hence 
providing a good indication that the periodic attractor 
remains essentially intact. 
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FIG. 4. Dependence of the average blowout variable (y) 
as a function of At, for fixed input signals. The parameters 
values are as in Fig. [[]. 

The above arguments and results demonstrate the 
qualitative differences between the responses of the on- 
off and the in-out dynamics to small input signals. In 
particular, the survival of the periodic orbit in the latter 
case acts to trap the incoming orbits and therefore blocks 
the possibility of supersensitivity in this case. We expect 
this picture to be common and thus supersensitivity to 
be absent in the generic non-skew product (in-out) set- 
tings. 
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FIG. 5. Dependence of the sensitivity S for a small in- 
put signal p as a function of At, together with the predicted 
scaling (^). The parameter values are as in Fig. |l|. 

To further substantiate this finding, we calculated the 
average blowout variable (y) in system (^J) for fixed input 
signals, as a function of At- The results are summarised 
in Fig. H, which show that for At > 0, both cases are 
relatively insensitive to input signals, whereas for At < 0, 
the on-off case is much more sensitive to input signals 
than the in-out case, with the latter dependence in very 



good agreement with our prediction (g). 

Finally we calculated the dependence of the sensitivity 
S for the in-out case as a function of At, with differ- 
ent input signals, as a function of At- The results are 
shown in Fig. pi together with our predicted expression 
(^), which show excellent agreement. 

To summarise, we have argued that the supersensitiv- 
ity found in j| for the case of skew product systems with 
on-off intermittency is absent in the more general set- 
ting of non-skew product systems, capable of displaying 
in-out intermittency. We have substantiated this claim 
both analytically and through extensive numerical simu- 
lations. We have also checked that the absence of super- 
sensitivity in the in-out case remains robust to changes 
in both the input signal (of the form psin(27rcc)) as well 
as to unbiased noise in the transverse direction (of the 
order of the input signal). 

The absence of supersensitivity for systems displaying 
in-out intermittency is important, particularly given that 
dynamical systems are generically expected to be of non- 
skew product type. 
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